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ABSTRACT 


Most gun propellant is stored at depots for a long time 
before it is used. While being stored, the quality of the gun 
propellant may deteriorate and become unstable. In an attempt 
to avoid disaster due to use of unstable gun propellant, 
accurate prediction of the safe shelf-life of gun propellant 
is necessary. The shelf-life estimation methods used 
currently for a group of similar gun propellant lots are 
based on a fixed effects regression model. This does not take 
into consideration the fact that samples from the same lot are 
more Similar than samples between lots. To capitalize on this 
lot-to-lot variation when estimating the shelf-life, first, a 
wanton: Sr sces regression model is developed. Secondly, a 
combined mixed effects model is estimated. The estimated 
model is then used to predict not only the shelf-life of a 
group of similar lots but also that of each individual lot of 
5"/54 NACO gun propellant stockpile. The results indicate 
that, first, the claimed shelf-life 1s not adequate and 
requires amendment. Next, the group shelf-life estimated can 
be relatively conservative compared to the individual shelf- 
lives. In view of potential opportunity loss due to safe 
individual lots being discarded, use of individual shelf-life 


is recommended. 
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EXECUTIVE SUMMARY 


Most gun propellant is stored at depots for a long time 
before it is used. While being stored, the quality of the gun 
propellant may deteriorate and become unstable. In an attempt 
to avoid disaster due to use of unstable gun propellant, 
accurate prediction of the safe shelf-life of gun propellant 
is necessary. Before 1914, real time storage inspection was 
used to estimate the safe shelf-life of gun propellant lots 
produced in the U.S.A. In an effort to reduce the length of 
time required for real time storage inspection, an Accelerated 
Aging Test (AAT) was introduced at Picatinny Arsenal in 1914. 
Under AAT, master samples taken from gun propellant lots are 
sequentially stored in heating chambers at 65.5 °C to measure 
the red nitrogen oxide fume time. Once the fume time is 
observed, it is related to the age of propellant stockpiles. 

The current methods to estimate the shelf-life of a group 
of similar gun propellants are based on a fixed effects linear 
regression model for the fume time. The shortcomings of the 
current estimation methods are such that the relationship 
between the fume time and propellant age may not necessarily 
be a linear one and the fixed effects model does not take into 
account the potential lot-to-lot variation. To consider 
potential lot-to-lot variation, a random effects regression 
model (RERM) is required. 

To understand the state-of-art estimation of the safe 
shelf-life of gun propellant, four technical reports are 
briefly reviewed. Furthermore, a preliminary analysis is 
conducted to understand the basic characteristics of the 
master samples and fleet return data sets of 5/54 NACO gun 
propellant provided by Naval Surface Warfare Center (NSWC). 
By analyzing the fitted models and validating the model 


Vil 


assumptions, the standard linear regression model based on the 
reduced data set is selected to further investigate a random 
effects regression model analysis. 

To consider a potential Jlot-to-lot variation when 
estimating the shelf-life, first the RERM is developed using 
a two-stage analysis. Next, a combined mixed effects model is 
estimated by employing both the maximum likelihood (ML) and 
restricted maximum likelihood (REML) methods. These 
estimators are used to define the shelf-life of the group of 
similar lots. The estimated model is used to predict not only 
the shelf-life of a group of the similar lots but also that 
for each individual lot by a shrinkage procedure. 

The results indicate that, first, the claimed shelf-life 
(35 years) for the NACO propellant lots is not adequate 
because it overestimates the shelf-life for NACO propellant 
lots by almost 10 years. Therefore, it requires amendment. 
Next, the group shelf-life estimated can be relatively 
conservative compared to the individual shelf-lives. In view 
of potential opportunity loss due to individual lots discarded 
when they are still safe, use of individual shelf-life is 
recommended over the group shelf-life when the management of 
individual lots is expected to cost less than that of bulk 
management. 

There is a nonlinear decreasing trend in fume time in some 
individual lots as the propellant ages. This raises questions 
about the approaches that use traditional linear least squares 
analysis of fume time data to predict safe shelf-life. This 
thesis recommends that a nonlinear regression analysis 


approach be conducted for some individual propellant lots. 
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I. INTRODUCTION 


A. BACKGROUND 

Most gun ammunition is stored at depots for a long time 
before being used. In order to avoid potential disaster due 
to storing unstable gun propellant, the shelf-life of a gun 
propellant stockpile should be accurately estimated. 
Typically, the shelf-life of gun propellant stockpile is 
defined as the age at which 5% of the gun propellant is 
unstable. The estimated shelf-life can also be used to 
efficiently minimize the ammunition management costs for 
distribution, maintenance, consumption and disposal. 

In the early 1900’s, real time storage inspection was 
used to estimate the safe shelf-life of gun propellant lots 
produced in the U.S.A. However, under the real time storage 
inspection, it took a long time to find visible changes in 
stability of propellent. In an effort to reduce the length of 
time required for the real time storage inspection, an 
Accelerated Aging Test (AAT) was introduced at Picatinny 
Arsenal in 1914 [Ref.3]. For AAT, master samples are 
collected for testing purposes from each gun propellant lot 
produced in the U.S.A. over the past sixty years. Each master 
sample consists of 5 pounds of propellant from each lot. 

Under the AAT, forty-five grams of propellant are taken 
from each master sample and are stored in heating chambers at 
65.5 °C until red nitrogen oxide fumes appear. The time it 
takes to observe such oxide fumes is called the fume time. 
This fume time is recorded with the corresponding propellant 
age. If the fume time falls below thirty days, the propellant 
is considered unsafe and the lot is recommended to be 
condemned. A fresh batch of propellant is then taken from 
ambient storage and the AAT is repeated. In general, it is 
assumed that fume time decreases as propellants age. In order 


to validate the results obtained from the master samples, the 











AAT is employed on the fleet-return data of gun propellant 
lots stored in the active depots. 

The classical methods used to estimate the shelf-life of 
a group of similar gun propellants is based on a fixed effects 
linear regression model. The fume time data are aggregated 
over a group of similar lots and a linear model for fume time 
is fitted against propellant age. The safe shelf-life of 
current methods is then estimated as the time period at which 
a 95% one sided lower prediction limit for the fume time curve 
intersects the acceptable lower specification level of thirty 
days of fume time. The shortcomings of the current estimation 
method are such that the relationship between the fume time 
and propellant age may not necessarily be a linear one and the 
fixed effects model does not take into account the potential 
lot-to-lot variation. Other relationships in addition to the 
linear model should be examined. In order to consider 
potential lot-to-lot variations, a random effects model is 
required. 

To understand the state-of-art related to the estimation 
of the safe shelf-life of gun propellant, four technical 


reports are briefly reviewed. 


B. LITERATURE REVIEW 

In the report entitled "Prediction of Safe Life of 
Propellants" [Ref. 1], the nature of 65.5°C surveillance test 
and the propellant chemical deterioration are briefly 
discussed. In this study, it is recognized that the 
measurement of residual stabilizer content offers the best 
means of establishing the stability potential of propellant. 
This report also discusses the results obtained from the 
standard artillery propellant when exposed to the aging test 
at various temperatures. | 

The measurements of residual stabilizer content versus 


time were used as the proxy for propellant deterioration. For 














interpretation of this deterioration phenomena, Berthelot’s 
Law was employed. It demonstrated that a family of straight 
lines can be plotted characterizing the length of time 
necessary at the various test temperatures to obtain a given 
variation of stabilizer content. Finally, by establishing 
realistic cut-off points regarding stabilizer content in a 
given storage temperature, propellant safe life was estimated. 

Secondly, the report entitled "Type Life Program History, 
Philosophy, Accomplishment, And Future" [Ref. 2] introduces 
the concept of the type life test. The type life test is 
derived from the AAT approach and takes into account a more 
realistic temperature. In this study, accelerated aging 
techniques are used to observe the aging mechanisms and 
characteristics of propulsion ordnance along with associated 
components. They are also use to determine how age affects 
the performance of the units. 

The type life temperature-time profile defined in this 
study is based on the "hot-month" concept, and differs from 
the AAT and its use of constant testing temperature. The 
concept is briefly stated as follows:(1)propellants are a 
chemical system; (2)chemical systems deteriorate more rapidly 
at elevated temperatures; and (3)in a normal year’s time, the 
deterioration occurs only during the summer months. 

Based on "hot-month" concept, the engineers establish a 
temperature-time profile termed the "compressed-ambient cycle" 
and consisting of 26 weeks. This profile was designed to 
Simulate the four seasons of the year and is related to the 
aging that normally occurs during magazine storage in the 
U.S.A. There 1s a 16-week period at 100 °F (38°C) representing 
a long severe summer time, two 3-week periods at 70 °F 
representing spring and autumn, and a 4-week period at 40 °F 
representing winter. The compressed-ambient cycle (26 weeks 
or 1/2 year) simulates one year of magazine storage; therefore 


there is a 2:1 (1.e.52 weeks:26 weeks) aging ratio. The main 


purpose of this report was to evaluate the effect of the loss 
of nitroglycerin from the propellant into the inhibitor. This 
loss results in propellant cracking and/or a reduction in 
propulsion performance due to a loss of energy. 

Thirdly, a report entitled "Statistical Analysis of 
ARDEC’S Fume Data Base" [Ref. 3] describes the procedure for 
AAT based on the master sample surveillance program in detail. 
Three types of propellants were used separately based on 
single and double base propellants. The single and double 
bases are terms of ammunition, distinguished by the methods of 
ignition. The single base propellant can be ignited by either 
mechanical or electrical means. The double base propellant 
can be ignited using both methods. In this study, fume times 
are plotted against sampling times for each propellant form. 
Not only is the first order linear relationship examined with 
this approach but also other relationships such as polynomial 
regressions, logarithmic and exponential transformation are 
considered to establish fume time behavior over the life of 
the propellant. 

Finally, a study entitled "Long-Term Stability of Navy 
Gun Propellants" [Ref. 4] was initiated to standardize the 
requirements of safety surveillance testing between the Army 
and Navy. This was necessary because current Navy 
surveillance testing includes a 65.5 °C oven fume test (AAT) 
on single-base propellants and double-base propellants with 
less than or equal to §<106 -‘Hatroglycerim content. 
Additionally, propellants must take more than 30 days to fume 
to be considered stable. In contrast, the Army’s disposition 
criteria is based on a stabilizer analysis of less than 0.2% 
nitroglycerin content. Because a substantial quantity of Navy 
propellant is being currently stored in Army depots, non- 
standardization complications commonly occur. This study was 
tasked by Naval Sea Systems Command(NSSC) to determine the 


critical level of propellant stabilizer for selected 











propellants of interest to the Navy and to the Marine Corps, 
to develop the best test methods to evaluate that level of 
criticality for Navy propellants, and to apply the data for 
standardization. 

The test was conducted on propellants removed from loaded 
ammunition. Each sample was subjected to chemical, kinetic, 
and physical properties analyses. As a result of these 
analyses, it was recommended to continue use of the 65.5 °C 
oven fume test because it is the procedure that provides the 
most complete coverage of propellant in storage worldwide. 
The Navy has a large historical database of fume data, and 
statistical analysis of the data provides an inexpensive, 
comprehensive means of insuring that Navy-developed propellant 
is safe for worldwide storage. 

In summary, the master sample surveillance test is 
recommended for continuous use to determine safety of the gun 
propellant in storage. It is proposed that the relationship 
between the fume time and propellant age may not be limited to 
a linear relationship. The management of data structure and 
the data transformation (e.g. logarithmic or reciprocal) is 
required. Additionally, no potential lot-to-lot variation in 
shelf-life estimation was examined in the four technical 


reports summarized here. 


C. THE SCOPE OF THESIS 

In this thesis, the master sample surveillance test data 
set is first examined using several fixed effects regression 
models in an attempt to find the relationship between the fume 
time and propellant age. Then, a random effects regression 
model is introduced to take into account a lot-to-lot 
variation in fitting the relationship between the fume time 
and propellant age. Finally, results are discussed and 


recommendations are made. 








II. FIXED EFFECTS REGRESSION MODEL 


A. DATA DESCRIPTION 

A preliminary analysis is conducted to understand the 
basic characteristics of the master sample data set provided 
by Naval Surface Warfare Center (NSWC). The data set was 
collected based on the oven test: data consisting of 29 lots of 
master samples of 5"/54 NACO gun propellant. They were 
manufactured at Badger Army Ammunition Plant (BAAP), Baraboo, 
Wisconsin between 1972 and 1974. The NACO is a generic term 
of NAvy COol representing a family of cool-burning, single- 
base propellants. These propellants are made from _ low- 
nitration nitrocellulose with a coolant, a stabilizer, a 
decoppering agent, and a flash suppressant added. 

First by subtracting the manufacturing date (see Appendix 
A) from the date when a 45g propellant batch was sent to the 
heating chamber, the age is computed for each propellant lot. 
The data of the fume time and propellant age is shown in 
computer program list in Appendix B. Both descriptive 
statistics of fume time and propellant age are displayed in 
Table 1. In summary, 221 observations of fume time and 
propellant age had been taken from 28 propellant lots. The 
maximum (1492 days) and minimum (19 days) of fume time are 
observed during 6813 days of propellant age. Next, in an 
attempt to find the relationship between the fume time and 
propellant age, the fume time is plotted against the 
propellant age in Figure C.1. It is observed that there is an 
unexpected increasing trend in fume time over propellant age 
up to 1100 days. This unexpected pattern may be due to the 
type of ammunition, shape of propellant container or caused by 


some unknown chemical reaction. [Ref. 3] 





Descriptive Statistics For Raw Data 


manufacture date 


number of repeated 


measurements 


222.¢L59 


fume time 


ropellant 
proP 2020.824 
age in days 





Table 1. Descriptive Statistics For Raw Data 


B. STANDARD LINEAR REGRESSION MODEL I 
Regardless of such unexpected pattern, according to the 
current practice, standard linear regression model I, is 


applied to fit the fume time (Y) against propellant age (X): 


Y=B, +B, X+e (1) 


where B, and B, represent the expected initial condition and 
deterioration rate of stability of gun propellant, and € is 


the random error which follows a Normal (0,77) distribution. 








The predicted value (¥Y,) at X=X. is 
VY, = B, + B, x, (2) 


where B, and B, are the ordinary least square estimators of By 


and B,. Subsequently a 90% prediction interval for a fume 


time at given propellant age (X.) is 
Y. +t (0.95 , N-2 ) SEPRED ( Y, | X,), (3) 


where N is the total number of observations, t(0.95,N-2) is 


the 95th percentile of the t distribution with N-2 degrees of 
freedom. SEPRED (¥Y, |X,) is the standard error of predicted Y, 


at X= Ky : 


SEPRED ( Y,|X,) =t [i++ (4) 


Based on this information, the group shelf-life (t,) can be 
estimated by checking the time period when the 95% one-sided 


lower prediction limit intersects 30 days of fume time: 
30 = B+ B, t,,- € (0.95, N-2) SEPRED ( ¥,|t,,). (5) 


Next, actual data obtained from the NACO propellant 
surveillance test is used to fit the model and to estimate the 
shelf-life. The fitted model based on (2) is given in Table 

















lack of £it test degrees of freedom: 200,19; 


















P value 






P value = 0.004 


not Normal 


Normal distribution test 









(Kolmogorov-Smirnov Test) 


estimated shelf-life 10400 days (28.49 years) 


Table 2. Fitted Standard Linear Regression Model I 






2 and displayed as a solid line in Figure C.1. A 90% two- 
sided prediction interval (solid curves) is overlaid in Figure 
C.1. The resultant shelf-life based on (5) turns out to be 
10400 days (28.49 years) as shown in Figure C.2. 


The residual analysis (r; = Y, - ¥;, for i=1...N) for 


standard linear regression model I is shown in Figure C.3. 
Residual pattern indicates nonlinearity and nonconstant 
variance. Also, the Kolmogorov-Smirnov test statistic shown 
in Table 2 indicates that the residuals do not follow a normal 
distribucion:. 

For further data analysis, the data set that covers from 
zero to 1100 days of propellant age is deleted, because the 


pattern is unstable during this early stage and this pattern 


EO 








is not expected to last longer than 1100 days. Furthermore, 
after 1100 days of propellant age, the fume time tends to 
decrease as propellants age and this pattern is more 
consistent. Therefore, the use of the remaining data set 
appears to be reasonable for the shelf-life estimation. The 
following three regression models are employed (1) standard 
linear regression model MII (delete observations that 
correspond to propellant age below 1100 days and use standard 
linear regression model I);(2) lJlinearizing transformation 
model I (transform propellant age to 1/ propellant age); and 
(3) linearizing transformation model II (transform fume time 
to log(fume time)). Graphical method precedes statistical 


analysis. 


C. STANDARD LINEAR REGRESSION MODEL II 

The standard linear regression model II is essentially 
the same as the standard linear regression model I(1). 
However, parameters are estimated based on a reduced data set. 
After deleting the early age, the fume time is plotted 
against propeliant age in Figure C.4, The results of model 
fit are shown in Table 3. The predictive value (middle solid 
jiné) and a 90% two-sided prediction anterval. for a.fume time 
(solid curves) are overlaid in Figure C.4. 

The residual analysis for the standard linear regression 
model II is shown in Figure C.5. In general, the plot of 
residuals against fitted values over fume time between zero 
day and 750 days (propellant age above 4500 days) does not 
show any systematic features and they appear to have a common 
variance. But, a minor nonrandom pattern is observed when the 
fume time exceeds 750 days (propellant age below 4500 days). 
Also, based on the p-value associated with Kolmogorov-Smirnov 
test shown in Table 3, we can conclude the residual of this 


model follows a normal distribution. 


a0 





Standard Linear Regression Model II 


lack of fit test degrees of freedom: 167,11; 
P value = 0.001 





? a 








Normal distribution test P value = 0.458 


(Kolmogorov-Smirnov Test) Normal 


estimated shelf-life 8469 days (23.20 years) 


Table 3. Fitted Standard Linear Regression Model II 


By checking the level of 95% one-sided lower prediction 
limit for a fume time that intersects 30 days of fume time in 
Figure C.6, the estimated group shelf-life turns out to be 


8469 days (23.20 years). 


D. LINEARIZING TRANSFORMATION MODEL I 
In linearizing transformation model I, it is assumed that 


the fume time is linearly related to reciprocal of propellant 


age: 


) + € (6) 


The fume time is plotted against propellant age in Figure C.7. 


12 

















The results of a model fit are shown in Table 4. The 


predicted value for the fume time ¥, = B, + B, ( = ) is shown 






Linearizing Transformation model I 
SE (B,) 17.420 


SE(B,) 51308.706 


lack of fit test 







degrees of freedom: 167,11; 







P value 






P value = 0.02948 


not Normal 


Normal distribution test 










(Kolmogorov-Smirnov Test) 


Table 4. Fitted Linearizing Transformation Model I 





as a middle curve and a 90% two-sided prediction interval for 
a fume time are overlaid in Figure C.7. 

It is interesting to note that the 95% one-sided lower 
prediction limit for predictive fume time converges to around 
400 days as the propellant age continuously increases. 
Therefore, the safe shelf-life of gun propellant based on 
linearizing transformation model I cannot be estimated or can 
be considered as infinity. 

The residual analysis for the linearizing transformation 


model I is shown in Figure C.8. In general, the plot of 


3 


residuals against fitted values over fume time between 400 
days and 800 days (propellant age above 4000 days) does not 
show any systematic features and they appear to have a common 
variance. But a nonrandom pattern is observed when the fume 
time exceeds 800 days (propellant age below 4000 days). This 
phenomenon is similar to the standard linear regression model 
II (see Figure C.5). Also, based on the p-value associated 
with Kolmogorov-Smirnov test shown in Table 4, we can conclude 
that the residual of this model does not follow a normal 


aq: stribution. 


E. LINEARIZING TRANSFORMATION MODEL II 
In linearizing transformation model II, it is assumed 


that the propellant age is linearly related to log fume time: 


LOG (Y) =B, +B, X+€é. (7) 


The fume time is plotted against propellant age in Figure C.9. 


The results of a model fit are shown in Table 5. The 
predicted value for the fume time YY, = EXP (B, + BLX,) anda 90% 


two-sided prediction interval for a fume time are overlaid in 
Figure C.9. 

The residual analysis for the linearizing transformation 
model II is shown in Figure C.10. In general, the plot of 
residuals against fitted values over fume time between Zero 
day and 6.6 days (propellant age above 4000 days) does not 
show any systematic features and they appear to have a common 
variance. But a nonrandom pattern is observed when the fume 
time exceeds 6.6 days (propellant age below 4000 days). This 
phenomenon is also similar to the previous two models (see 
Figure C.S and Figure C.8). Also, based on the p-value 


associated with Kolmogorov-Smirnov test shown in Table 5, we 


14 











Linearizing Transformation model II 


lack of fit test © degrees of freedom: 167,11; 


P value 


Normal distribution test P value = 1.0E-11 


(Kolmogorov-Smirnov Test) not Normal 


estimated shelf-life 17400 days (47.67 years) 


Table 5. Fitted Linearizing Transformation Model II 





can conclude that the residual of this model does not follow 
a normal distribution. 

By checking the level of the 95% one-sided lower 
prediction limit of a fume time that intersects 30 days of 
fume time in Figure C.11, the estimated group shelf-life turns 
out to be 17400 days (47.67 years). | 

In summary, by analyzing the fitted models’ and 
validating the model assumptions of the various fixed effects 
regression models, the effective distribution data model 
appears to be the standard linear regression model II based on 
the reduced data. This model is compared to the model fit 


based on the fleet return data. 


i 





F. STANDARD REGRESSION MODEL BASED ON FLEET RETURN DATA 
This model is based on the oven test data consisting of 

28 lots of fleet return data of 5"/54 NACO gun propellant. 

The descriptive statistics of fume time and propellant age are 


displayed in Table 6. In summary, 204 observations of fume 


manufacture date 


number of repeated 


measurements 


fume time 


in days 


propellant 


age in days 





Table 6. Descriptive Statistics For Fleet Return Data 


time and propellant age had been taken from 28 propellant 
lots. The maximum (993 days) and minimum (479 days) of fume 
time are observed during 5074.5 days of propellant age. The 
fume time is plotted against propellant age in Figure C.12. 


The results of model fit which is based on standard linear 
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degrees of freedom: 127,73; 
P value = 0.0000 









Normal distribution test P value = 0.0040 






(Kolmogorov-Smirnov) not Normal 


estimated shelf-life 15700 days (43 years) 


Table 7. Fitted standard Linear Regression Model 
Based on Fleet Return Data 


regression model I are shown in Table 7. The predicted value 
(middle solid line) and a 90% prediction interval for a fume 
time (solid curves) are overlaid in Figure C.12. The 
resultant shelf-life based on (5) of this model turns out to 
be 15700 days (43 years) as shown in Figure C.13. 

The residual analysis for this model is shown in Figure 
C.14. There is a right-opening megaphone pattern in Figure 
C.14. This residual pattern observed shows nonconstant 
variance. Also, the Kolmogorov-Smirnov test statistic shown 
in Table 7, indicates that the residuals of this model do not 
follow a normal distribution. Two remedies for nonconstant 
variance (Weisberg [Ref. 8]) are made using weighted least 


square estimation and power transformation: 


Le 











1. Weighted Least Square Model 

The first remedy is the weighted least square model. The 
variance of the random error in the standard linear regression 
model I (1) appears to increase as X increases, that is, 7T=KX 


where K is a positive constant. Therefore, the model (1) is 


modified as follows: 


(8) 


The residual analysis for this model is shown in Figure C.15. 
It still shows a right-opening megaphone pattern in Figure 
C.15 and the suggested weighted least square method is not 
adequate. 

2. Power Transformation Model 

The second remedy is to transform the response Y via a 
variance stabilizing transformation. By checking the residual 
pattern of several transformation models, the following 
transformation model appears to satisfy the assumption of 


common variance (see Figure C.16): 


ye = 3. ee A ee) (9) 


The fume time is plotted against propellant age in Figure 
C.17. The results of a model fit are shown in Table 8. The 
= a - -2 
prediction fume time Y= (B, + B, X,) ° and a 90% prediction 

interval for a fume time are overlaid in Figure C.17. 
The p-value associated with Kolmogorov-Smirnov test of 
residual of model (9) indicates that the residual of this 


model follows a normal distribution. 
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Lack -0f E21 -test degrees of freedom: 127,73; 
P value = 0.001 

















Normal distribution test P value = 0.18964 





(Kolmogorov-Smirnov) Normal 


estimated shelf-life 


Table 8. Fitted Power Transformation Model 


The shelf-life is estimated based on the following 


fFOLrMii La: 


-1 (10) 
30 2 [By By bj 6 (- 0595: 5. N-2) SEPRED( Y,.| tay 1. “a 


The estimated group shelf-life turns out to be saben Gast wo 
This shelf-life estimation based on fixed effects 
regression models does not take into account the lot-to-lot 
variation. Apparent ly potential lot-to-lot variation is 
supported by individually fitted models based on the Standard 


Linear Regression Model I, II and fleet return data as shown 


9 











in Figures C.18, C.19, C.20. Additionally, this variation can 
be seen by checking linearizing transformation model I after 
smoothing the data following Lowess procedure [Ref. 9] as 


shown in Figure C.21. 
In the next chapter, standard linear regression model II 


will be further investigated for a random effects regression 


model analysis. 
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TII. RANDOM EFFECTS REGRESSION MODEL 


A. INTRODUCTION 

The current methods used to estimate the shelf-life of 
gun propellant is based on the standard linear regression 
model I. In this fixed effects linear regression model, the 
expected initial fume time (B,) and deterioration rate (B,) are 
assumed to be constant for all gun propellant lots. However, 
in practice, the performance varies among lots even though a 
group of lots is manufactured based on the same technique. 
This kind of phenomenon can be observed in the individually 
fitted. models in Piqure. €.168 to .C.21. The lot-to-lot 
variation is observed not only in the intercepts (expected 
initial condition) but also in slopes (deterioration rate) 
among a group of similar lots. The current methods does not 
take into account this kind of constant. In order to consider 
potential lot-to-lot variation within the group, a random 
effects (RERM) (or mixed effects) regression model is 
necessary. 

In this chapter, first, the RERM will be introduced based 
on a two-stage analysis. Next, a combined two-stage model 
(mixed effects model) will be fitted based on the reduced data 
of the master samples provided by NSWC. In order to obtain 
the estimators of unknown parameters in mixed effects model, 
the maximum likelihood estimation (MLE) and restricted maximum 
likelihood (REML) will be employed. These estimators are, 
then, used to defined the shelf-life of the group of similar 
lots that takes into account the lot-to-lot variation. 
Furthermore, a shrinkage procedure will also be introduced to 


estimate the shelf-life of an individual lot. 


C70 











B. MODEL 
1. Within-Lot Model 


In order to estimate the deteriorating patterns of gun 


propellant lots over storage time, each 45g batch is taken 
from 5 pound master sample of lot i and its ageing is 
accelerated at the heating chamber at 65.5 °C. The fume time 
(yj) of lot i (i=1...N) is recorded with the corresponding 
propellant age (t,) (j=1...n,;) whenever the red nitrogen oxide 
fume appears. Once repeated measurements are obtained, we 
attempt to find the relationship between the fume time and 
propellant age. The following within-lot model is used to 
describe the relationship between the fume time and the age of 
a lot: 
For i-1...N and J=H1l...n,; 


Viz = Moy + Oy; Cy, + €y; (11) 


where a), and a, represents the expected initial condition and 
deterioration rate of the stability of propellant lot i, and 
€é; is the random error which 1s assumed to follow a 
independent N(0, 7’). 

2. Between-Lot Model 

In order to take into account the lot-to-lot variation 
among intercepts (expected initial condition) and slopes 
(deterioration rate) of a group of lots, we assume the 


following between-lot model: 


Gos = Wy Yo + So, 
(12) 
Os: = W; Yi 3 4; 


Zz 











where w, is a 1 x k vector representing values of K factors 
that may cause varying a, and a; . The yo and y; are k x 1 


vectors of regression coefficients for 4d, 


1 


and Qj, 
respectively. The 6) and 6,, are assumed to follow independent 


2 2 
Oo Oo1 


N(0,)), where ee Note that for AAT, the master 


2 2 
Oo1 01 


samples are put in the constant temperature condition (65.5 
°C) and w, in this case is 1. 

When a) and a, in (11) are substituted with those in 
(12), a combined mixed effects model and is formulated as 
follows: 


POY D=i,..Noand. jelevs.n, 


= ( Wy Yo + We Ya Cay ) + ( Oo; + 84; Ci ) 4 ©ij (13) 


o 


z a i 
(mes wy tig) (Te # Ca, ta) (Be ] eu 


This mixed effects model can be re-written using matrix 


notation: 


y=xXPB + Z2vtre, (14) 


where y is a n,;x1 vector of yj; 6 is a 2Kx1 coefficient vector(y, , Va) 


N 


associated with a >) 1;x2K matrix of X consisting of 
on 


Ze 











wt v is a 2njx1 vector of (5); , 5,,)’ associated with 


(W, , W, iz) ; 


N 
a known )°n,x2n, block diagonal matrix Z of (1, €;;), and e 
L?t 


N 
is an >, 23x1 vector of random error. Assume that uv and e 
L=1 


are uncorrelated and have expections 0 and variances G and R, 
respectively, where both G and R are nonsingular. applying 


these assumptions to (14) leads to 


(15) 


When G and R (and hence V) are known, estimates of both 


6 and the realized value of v are 


B= [% + 91]! = (x’V7x) *x'v ty 
(16) 


O=[0,,... Oy] =[8,,,815, - + Soy Syl =GZ/V > (y-XB) 


respectively, where 6 is the best linear unbiased estimator 


(BLUE) of 8, and % is the best linear unbiased predictor 
(BLUP) of ov. They can also be obtained from Henderson’s 
(1984) mixed model equation [Ref. 10]: 
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X’/R4X eZ | | 6 | (17) 


Z'R ix Z'’R-1Z+G71 








where B = [Yo 5 47 and ¢ = [So BOs aa Ones Sia When G and R 


(and hence V) are unknown, the estimator of V can be obtained 


by maximizing the logarithm of likelihood function [Ref. 10] 


log L=- N logan -=1og|v|- (y-XB)! V-t (y-xB) . (18) 


The ML estimators V and ff can be obtained by solving the 


following equations simultaneously 


ROW Sp eae ye ey (19) 
and 
S—1 4 [orl lft ; (20) 
ELV eZ) = AVR AB) IV ZA) 0 Or 7 =05 2 jaa, 
N 
where Z is a n;x2n,; block diagonal matrix of (1,t,) and 
Lei. 
N N 


priv +22.) a6 Gvace. of the » 2x), 2; matrix (¥72Z,Z:) that 
1=1 1=1 


denotes the sum of the diagonal entries. An algebraically 


Simpler expression for (20) is derived by defining 


Z5 














P=v - wix(x'/v1x) 1x’v-, (21) 


Then from (19) it is clear that for P being P with V replaced 


ae 


by 


Oty - xB) = By, (22) 


so that the ML estimators of V and 6 also can be re-written 


by solving the following ML equations simultaneously 


SOE SOR ae Pe (23) 
and 
te ORV Zi) Te Re PLE Ft (24) 


where | tr (072,23) }_, is (N+1)x1l column vector with elements 


which are the trace of the matrix (V2.7). FOr 204.5% oN, 


respectively. 

Also, the estimator of V can be obtained by restricted 
maximum likelihood (REML). A basic idea of REML estimation is 
that of estimating variance components based on residuals 
calculated after fitting by ordinary least squares just the 
fixed effects part of the model. The REML estimator can also 
by viewed as maximizing a marginal likelihood. The REML 


equations can therefore be derived from the ML equations of 


(24). By making suitable replacements, the REML estimator V 
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can be obtained by solving the following equation: 


{ tr [ (K/VK) *K'Z2'K l}i.5 = 
(25) 
{ y/K(K/OK) *K'ZZ'K( K'OK) *K'y }F, 


with P = K(K’VK)'K’ where K’X=0. 
When the resultant V is found by solving the ML equations 


((19),(20)) or REML equation (25), the Henderson’s mixed model 


equation can be re-written as follows: 


X’/R 1x Rg 


: a (26) 
VAS ie 4 Zip + Z+G 





OF 


Finally, B =[f. QJ and 9 = [8),,8,,,...89y,8,,) substitutes 








B = nee aia and v = [S17 9i1, ve Son San) in (14) as follows: 


P=xXP+29. (27) 


In fume time analysis, the test temperature condition is 
constant and W, in this analysis is 1. Therefore, 


X= (W., Wt) = (1 1 t;) and 4i= (1 t Cj) e Then (27) becomes 


Pe | 


xB + Z0= ZB + ZO 


> 
T 


= ZB + Z(GZ/V" (y-ZB) ) 


28 +Z6Z'Vy-ZGZ'V"* zB (28) 


(1-262'7 7) ZB +(ZEz'V) y 


(I-H) B(y) +Hy 


where H=Z6Z'/7+ is the weighted of mean of Ely] and y. The 
resulting Y is called a shrinkage estimator. 


The predicted value of fume time (¥,;;) of lot i at tj can 


be estimated as 


A 90% prediction interval for a fume time at given propellant 


age ts 2s 


¥i;, +t (0.95 , N-2 ) SEPRED ( Hz, | C35) (30) 


where 


1 
SEPRED () s(G)a( (ive) Vib) Grte) ee) VO) Git) (7; 
| (31) 


and V (B,;) and V (9,) is the i" block diagonal matrices of 


28 














(32) 


In order to estimate deteriorating patterns of a group of 
Similar lots, one would be interested in the shelf-life of the 
Sroun- of cimilar-Lots: (tq that cakes ante account the: 1oc-t6- 
lot variation. It can be obtained as follows by only 


considering the group effects: 


30°] 209 P= 6 0 0495: 4. a2 ) SEPRED ee. Yh Ea3.) 


(33) 
0 


= (1, ten( 9°] -£ (0.95,N-2) SEPRED group) (¥ | ty2) 


where SEPRED,roup) ( ¥ | €3;) is the standard error of predicted 


y of a group given by 


1 
SEPRED ook len) =) (ipo) VR (ist) aes 8 


The group shelf-life can be compared to the individual 
shelf-life (t,,) based on the shrinkage estimation (28). It 
is defined as the time period at which the expected fume time 


curve intersects the 30 days of fume time: 


ad 








307= xB + Z 
(35) 


=( I - Z6zZ'0-" ) ZB + zG2Zz'F ty 


where I is a nxl vector of 1’s and Z is a n,;x1 vector of 


(1, Cea) 


C. DATA ANALYSIS 

Standard linear regression model II based on the reduced 
data is further investigated for the random effects regression 
model analysis. As explained earlier, under the AAT, the fume 
time (y;) of each lot is repeatedly measured five to nine 
times in an accelerated condition (65.5 °C) during the last 20 
years. The age (t;) for each of the lot is computed by 
subtracting the manufacturing date (see Appendix A) from the 
date in which when a 45g propellant batch was sent to the 
heating chamber. 

First, in order to estimate the safe shelf-life of a 
group of similar gun propellants based on the random effects 
regression model, the within-lot and between-lot models were 
combined to be a mixed effects model (14). Since the G and R 
(hence V) are unknown (15), the estimators of G and R can be 


obtained by a restricted maximum likelihood (REML) (25): 


ee Go Go, 0.07079100 -0.00450653 
GO MAE) | ae “aol I at OORSOESS:. 20 
62, 6? 3 0.00030661 
and R = Var(e) = 0.05926819. Unknown parameters 8 and v can 
in turn be estimated by Henderson’s mixed equation (26) : (1) 
B= [%,91]' = [3.81021477, -0.14105394] where se[f,]= 
[0.07176290] and se[f,]=[0.00534921]; (2) 0 = [89,,84,,-- Soy, bin) 


and the corresponding standard errors are shown in Table 9. 
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AGE LOTID 11192 -~Q0.01170315 0.01302592 
INTERCEPT HOTLD:- LLEes 0.19615700 O~18504953 
AGE LOLLTD: 1193 =0, 01161353 7 0.01297891 
INTERCEPT LOTID 11194 0.776200 0.18452648 
AGE LOTID 11194 -0.01152574 0.01310043 
INTERCEPT LOTID: 11195 0.16493381 0.18380045 

LOTID. 11295 ~0.01078221 O, 01291217 


AGE 


Table 9. % And Std Error of 3 


Subsequently, the safe shelf-life of the group of similar lots 


(ta 
obtained as follows: 


that takes into account the lot-to-lot variation can be 


3.81021477 : 36) 
30° (1, te) ( ¢ Sarosgog) (0+: 95¢ 8-2) SEPRED group (7 | te! 


Figure C.22 shows the resulting shelf-life is 8405.12 days 
(23.02 years). 

This result can be compared to the one obtain by the 
(5), which is 8469 days Note 


that the current method without taking into account the 


current method (23.20 years). 
potential lot-to-lot variation provides a liberal shelf-life 
that is 64 days longer than the random effects regression 
model. 

Finally, individual shelf-life of lot i (t,,) based on the 
(28) 
period at which the expected fume time curve intersects the 30 
days of fume time (see Table 10). 
shelf-life (t,,) varies from 27.998 to 24.942 years. 


shrinkage estimation is computed by checking the time 
The resulting individual 
They all 
exceed the group shelf-life based on either the fixed effects 
(23.20 years) or the random effects (23.02 years) regression 
model . 

The group shelf-lives based on the prediction limit for 
the fume time can be compared to those based on the confidence 


limit for the expected value. They turn out to be 25.66 and 
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OBS LOTID shelf-life shelf-life shelf-life 
(shrinkage (confidence (prediction 


estimation) limit) Limit) 
1 11142 26.914 24.520 23a <4t27 
2 LLLas 21576 42o40i31 235974 
2. 17144 23;.056 25.411 24.315 
4 11145 26-691 24365 2229 
S “LLlas 26. 951 24.726 23.424 
6 11147 Zeit 24.729 235561. 
7 sAdITAasS 27.292 24.863 23;698 
8 11149 25.654 23%:630 222859 
& Alviso 272167 24.726 23.2561 
RO. Sata 27685 24.794 23.4629 
is. 22256 27.074 24.589 23.356 
LA Lis 5'7 242998 25.411 24.246 
L323 ,libos 26.433 24.178 7 ye a 0 
L4- Li259 A445. A ae 234422 
15: -do60 27.600 25/3069 25.6904 
16. List 262715 Zo. 308 Deon Ae 
Ly sLLee2 Z72576 25.068 23.4972 
18 11183 254563 2531 OS 22-4662 
19 11184. 255534 23629 225502 
20° -LLITS5 24.942 25 ode 225194. 
Zi. Lites 25.526 23.4024 ZEaxD IS 
22. ALALeg 25.603 23 1690 22.2624 
22 AA90 25:22:75 23.493 22.459 
24 11191 Vio ane pots. 2526.32 222.578 
Z2>:. LL 2 Zoe taZ 22.055 ZA hoo 
26. 11193 2536: /2 234771 225658 
27 11194 Zo. 590 234561 22.5500 
BO Le Zo G39 222098 Dae oN!” 
Table 10. Individual Shelf-Life of Propellant Lot 


24.314 years, which are longer than their counterparts based 


on the prediction limit. 


The individual shelf-lives can also be calculated by 


checking the time period when 95% one-sided lower prediction 


or confidence limits intersect 30 days of fume time. The 


results are shown in Table 10. 


About 42% of individual shelf- 


lives based on the prediction limit turned out to be shorter 
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than the group shelf-life (23.02 years). The rest are as same 
as the shrinkage estimation that all exceeds the group shelf- 
life based on either the fixed and random effects regression 
model. 

All calculations are done by using PROC MIXED and PROC 
REG of a statistical package SAS [Ref.12]. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Most gun propellant is stored at depots for a long time 
before it is used. In storage, the quality of the gun 
propellant may deteriorate and become unstable. In an attempt 
to avoid disaster due to use of unstable gun propellant, 
accurate prediction of the safe shelf-life of gun propellant 
is necessary. The shelf-life estimation methods used 
currently for a group of similar gun propellant lots are based 
on a fixed effects regression model. This does not take into 
consideration the fact that samples from the same lot are more 
Similar than samples between lots. To capitalize on this lot- 
to-lot variation when estimating the shelf-life, first, a 
random effects regression model (RERM) is developed. 
Secondly, a combined mixed effects model is estimated. The 
estimated model is then used to predict not only the shelf- 
life of a group of similar lots but also that of each 
individual lot of 5"/54 NACO gun propellant stockpile. 
Finally, results of the data analysis are discussed and 
recommendations are made. 

The claimed shelf-life of the NACO propellant lots 
provided by NSWC is 35 years. It 1s interesting to note that 
when the confidence limit method is used instead of the 
prediction limit method, the standard linear regression model 
I provides 35 years as the estimated shelf-life. However, the 
assumptions for standard linear regression model I were not 
validated. Therefore, the use of the claimed shelf-life (35 
years) is not supported and needs to be amended. 

The group shelf-life estimated based on a random effects 
regression model turns out to be 23.02 years. A fixed effects 
regression model based on the standard linear regression model 
II provides a slightly longer shelf-life (23.20 years) than 
23.02 years. The individual shelf-lives based on the 
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shrinkage estimation vary from 27.998 to 24.942 years which 
all exceed the group shelf-lives obtained from both the fixed 
and the random effects regression models. It indicates that 
group shelf-life estimated based on the fixed or the random 
effects regression models can be relatively conservative 
compared to the individual shelf-lives. When the group shelf- 
life (23.02 years) is applied to the maintenance of unstable 
stockpiles, there is a chance that all individual lots would 
be discarded even though they are still safe. In view of 
such opportunity loss, use of individual shelf-life is 
recommended over the group shelf-life when the management of 


individual lot is expected to cost less than that of the bulk 


management. 


B. RECOMMENDATIONS 

In this thesis, a particular propellant stockpile data 
set is analyzed. Although the approaches proposed can be 
applied to any other types of propellant, the model selection 
may be dictated by the nature of data set. 

There is a nonlinear decreasing trend in fume time in 
some individual lots as the propellant ages. This raises 
questions about the approaches that use traditional linear 
least squares analysis of fume time data to predict safe 
shelf-life. This thesis recommends that the nonlinear 


regression analysis approach be conducted for some individual 


propellant lots. 
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The following computer program is a SAS code. 


OPTIONS LS=80; 


DATA ALL1; 
INPUT OB 
ID= N ; 
CARDS; 
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APPENDIX 
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PROC 
PROC 


DATA 
MERGE 
IF ID 
IF ID 
PROC 
PROC 
PROC 


PROC 


PROC 


167 26. LiLigs 14.1219 702 
£88. 26 ‘-L1LL93 74.1219 703 
Lao. -26> “12193 74.1219 710 
LOO. 26: “L1bss 74.1219 726 
Lol. 26 D2EL9s 74.1219 918 
192 26 11193 12.1219 Tid 
£93 26: ALLS 74.1219 f gPRS 
194 26 11193 74.1219 726 
£95. 26. Lites 74.1219 731 
196 $26 12193 74.1219 697 
EO? 26 <Lit9s 1441219 700 
£98 26 -lt1935 74.1219 716 
£99.26: -<Lilgs 74.1219 654 
200° 26 1119s 74.1219 646 
201 27 11194 P4sL2i9 647 
202 27 11194 74.1219 688 
203 -26 ALT9S T4c1219 622 
204 28 11195 fe ee ee 609 


SORT; BY INDEX; 
PRINT DATA= ALL1; VAR 


/* 
NEW ; 
ALLI ALLZ2 ; BY INDEX; 

= LAG(ID) THEN EAGE = .; 


INDEX EAGE EDTF EYAGE EYDTF; 


COMBINE MASTER SAMPLE AND 


= LAG(ID) THEN EDTF = .; 
PRINT DATA=NEW; VAR INDEX MDATE DTF AGE EDTF EAGE; 
UNIVARIATE; 
PLOT; 
PLOT DTF*AGE='A’' EDTF*EAGE=’ *!’ ; 
PLOT; 
PLOT DTF*AGE='A’ EDTF*EAGE=' *’ /OVERLAY; 
/* STANDARD LINEAR REGRESSION MODEL I 
REG ; MODEL DTF = AGE ; 
/* STANDARD LINEAR REGRESSION MODEL II 
/* ( DELET EAGE BELOW 1100 DAYS ) 


DATA ALL3(DROP= OB ID EYAGE EYDTF);SET ALL1; 


DATA ALL4;SET ALL3; 
IF EAGE LE 1100 THEN DELETE; 


PROC 


PROC 


OUTPUT OUT=OU0T1 P= 


PROC 


PLOT ; 

PLOT EDTF*EAGE=% *' ; ? 
REG ; MODEL EDTF= EAGE; 
PRED R=RES; 
PRINT; 


DATA ALL5 (DROP= RES ) ;MERGE ALL4 OUT1; BY INDEX; 


PROC 


DATA ALL6 (DROP=OBS) ; 


PRINT; 


SET ALL5;BY INDEX; 


W= (EAGE-4146.15) * (EAGE-4146.15) ; 


SXX= 


SEPRED = 


2521553* (180-1) ; 
98.295 * SORT(1+(1/180) +(W/SXX) ) ; 


L95=PRED-1.645*SEPRED; 
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U95=PRED+1.645*SEPRED; 


EAGE30=30; 
PROC PRINT; 
PROC PLOT ; 
PLOT EDTF*EAGE='A’ PRED*EAGE=’-' L95*EAGE=’.’ 
U95*EAGE='.’ EAGE30*EAGE=' -' /OVERLAY; 


PLOT PRED*EAGE='’*' L9O5*EAGE=’.’ 
U95*EAGE='’.'’ EAGE30*EAGE=' -’ /OVERLAY; 
DATA EXTEN] ; 
DO DAGE = 8000 TO 12000 BY 100; 
PRED=1045.510-0.07*EAGE; 


W= (EAGE-4146.15) * (EAGE-4146.15) ; 

SXX= 2521553* (180-1); 

SEPRED = 98.295 * SQORT(1+(1/180) + (W/SXX) ) ; 
L95=PRED-1.645*SEPRED; 
U95=PRED+1.645*SEPRED; 


OUTPUT; 

END ; 

PROC PRINT DATA=EXTEN1;VAR EAGE PRED L95 U95; 
/* LINEARIZING TRANSFORMATION MODEL I * / 
/* (TRANSFER EAGE TO 1/EAGE) * / 


DATA ALL7 ; SET ALL1; 
IF EAGE LE 1100 THEN DELETE; 


Y=1/EAGE; 
PROC REG ; MODEL EDTF= Y; 
OUTPUT OUT=OUT2 P= PRED R=RES STDI=SEP; 


PROC PLOT; 
PLOT EDTF*EAGE; 
DATA ALL8 (DROP= RES ) ;MERGE ALL7 OUT2; BY INDEX; 


PROC UNIVARIATE; 
DATA ALLY (DROP=OBS); SET ALL8;BY INDEX; 


W= (Y-0.000297) * (Y-0.000297) ; 

SXX= 0.000164*0.000164* (180-1) ; 

SEPRED = 112.854* SORT(1+(1/180) + (W/SXX) ) ; 
L95=PRED-1.645*SEPRED; 
U95=PRED+1.645*SEPRED; 
PROC PRINT DATA=ALL9 ;VAR INDEX EAGE EDTF PRED U95 L95; 
PROC PLOT ; 

PLOT EDTF*Y=’A’ PRED*EAGE='-'’ L95*EAGE=’ .' 

USS*EAGE='.’ /OVERLAY; 


PLOT PRED*EAGE='*’ LO5*EAGE=’ .’ 
U95*EAGE='.’ /OVERLAY; 


DATA EXTENZ2; 
DO EAGE=15000 TO 35000 BY 100; 
PRED=411.7374+1315827* (1/EAGE) ; 


W= ((1/EAGE) -0.000297) * ((1/EAGE) -0.000297) ; 
SXX= 0.000164*0.000164* (180-1) ; 

SEPRED = 112.854* SQRT(1+(1/180) +(W/SXxX) ) ; 
L95=PRED-1.645*SEPRED; 
U95=PRED+1.645*SHPRED; 


4°] 











OUTPUT ; 


END ; 
PROC PRINT DATA=EXTEN2; VAR EAGE PRED L95 U95; 


[* LINEARIZING TRANSFORMATION MODEL II 
/* (TRANSFER EDTF TO LOG(EDTF) ) 


DATA ALL10; SET ALL1; 
IF EAGE LE 1100 THEN DELETE; 


EDTFLOG = LOG(EDTF); 
PROC REG ; MODEL EDTFLOG = EAGE; 
OUTPUT OUT=OUT3 P= PRED R=RES STDI=SEP; 


PROC PLOT; 
PLOT EDTF*EAGE ; 
DATA ALL11(DROP= RES ) ;MERGE ALL10 OUT3; BY INDEX; 


PROC UNIVARIATE; 
DATA ALL12(DROP=OBS); SET ALL11;BY INDEX; 


W= (EAGE-4146) * (EAGE-4146) ; 

SXX= 2521553* (180-1) ; 

SEPRED = 0.2742 * SQRT(1+(1/180) +(W/SXX)); 
L95=PRED-1.645*SEPRED; 
U95=PRED+1.645*SEPRED ; 


PROC PRINT DATA=ALL12 ;VAR INDEX EAGE EDTF PRED L95 U95; 


PROC PLOT ; 
PLOT EDTFLOG*EAGE='A’ PRED*EAGE=’-' L95*EAGE='.’ 


U95*EAGE='.’ /OVERLAY; 


PLOT PRED*EAGE='*' L9O5*EAGE=‘' .’ 
U95*EAGE=’.’ /OVERLAY; 


DATA EXTEN3; 
DO EAGE=14000 TO 20000 BY 100; 
EDTFH=EXP (7.3787-0.000197*EAGE) ; 
PRED =7.3787-0.000197*EAGE; 


W= (EAGE-4146) * (EAGE-4146) ; 

SXX= 2521553* (180-1); 

SEPRED = 0.2742 * SORT(1+(1/180) +(W/SXX) ) ; 
L95=PRED-1.645*SEPRED; 

U95=PRED+1.645*SEPRED; 

EXPL9S5=EXP (PRED-1.645*SEPRED) ; 

EXPU95=EXP (PRED+1.645*SEPRED) ; 

OUTPUT ; 

END; 

PROC PRINT DATA=EXTEN3;VAR EAGE PRED L95 U95 EXPL95; 


/* FLEET RETURN DATA * / 


DATA ALL14 ;SET ALL2; 
IF DTF GE 1200 THEN DELETE; 


PROC PLOT; 
PuUOT DIP*CAGE={"*0 3 
PROC REG ; MODEL DTF= AGE; 
OUTPUT OUT=OUTS5 P= PRED R=RES STDI=SEP; 
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PROC PRINT; 

PROC UNIVARIATE; 

DATA ALL15(DROP= RES ) ;MERGE ALL14 OUTS; BY INDEX; 
PROC PRINT DATA = ALL15; VAR INDEX DTF AGE PRED; 
PROC PRINT ; 


DATA ALL16; SET ALL15;BY INDEX; 
W= (AGE-2821.968) * (AGE-2821.968) ; 


SXX= 1928683* (202-1); 
SEPRED = 61.32 * SORT(1+(1/202)+(W/SXX) ) ; 
L95=PRED-1.645*SEPRED; 
U95=PRED+1.645*SEPRED; 


AGE30=30; 
PROC PRINT; VAR INDEX AGE DTF PRED U95 L95; 
PROC PLOT ; 

PLOT DTF*AGE='A’ PRED*AGE='-' L95*AGE=’.’ 


U95*AGE='.'’ AGE30*AGE=’ -’ /OVERLAY; 


PLOT PRED*AGE='*' L9O5*AGE=’ .’ 
U95*AGE='.' AGE30*AGE=’ -’ /OVERLAY; 


DATA EXTENS; 
DO AGE=10000 TO 16400 BY 100; 
PRED= 823.375-0.0429 *AGE; 


W= (AGE-2821.968) * (AGE-2821.968); 

SXX= 1928683* (202-1); 

SEPRED = 61.32 * SQRT(14+(1/202) +(W/SXX)); 
L95=PRED-1.645*SEPRED; 
TEST=PRED-1.699*SEPRED; 
U95=PRED+1.645*SEPRED; 

AGE30=30; 


OUTPUT ; 
END; 


PROC PRINT DATA= EXTENS5;VAR AGE PRED L95 U95 TEST; 
PROC PLOT; 


PLOT L95*AGE=! *’ AGE30*AGE=' .’ /OVERLAY; 
/* WEIGHTED LEAST SQUARE MODEL 
/* (DTF/AGE VS 1/AGE) 


DATA ALL20 ;SET ALL2; 
IF DTF GE 1200 THEN DELETE; 
NEWDTF=DTF/AGE ; 
NEWAGE=1/AGE; 
PROC REG ; MODEL NEWDTF= NEWAGE; 
OUTPUT OUT=O0T2Z1 P= PRED R=RES STDI=SEP ; 
PROC PRINT; 
PROC UNIVARIATE; 
DATA ALL22;MERGE ALL20 OUT21;BY INDEX; 
U9S. = PRED+1.645*SEP; 
L95 = PRED-1.645*SEP; 
PROC PRINT DATA=ALL22;VAR AGE DTF PRED L95 U95 SEP; 


PROC PLOT; 
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PLOT RES* PREO='*! ; 
PROC PLOT; 
PLOT L9S5 *AGE='*' PRED*AGE='’*’ U95*AGE=' *’ /OVERLAY; 


/* POWER TRANSFORMATION MODEL * / 
/* (DTF**-3 VS AGE) * / 
DATA ALL23 ;SET ALL2; 
IF DTF GE 1200 THEN DELETE; 
NEWDTF=DTF** (-3) ; 
PROC REG ; MODEL NEWDTF= AGE; 
OUTPUT OUT=OUT24 P= PRED R=RES STDI=SEP ; 
PROC PRINT; 
PROC UNIVARIATE; 
DATA ALL25;MERGE ALL23 OUT24;BY INDEX; 
U95 = PRED+1.645*SEP; 
L95 = PRED-1.645*SEP; 
PROC PRINT DATA=ALL25;VAR AGE DTF PRED L95 U95 SEP; 
PROC PLOT; 
PLOT RES* PRED=’*’' ; 
PROC PLOT; 
PLOT L95 *AGE='*’ PRED*AGE='*’ U95*AGE='*' /OVERLAY; 


DATA FIM4EXT; 
DO AGE=900000 TO 1000000 BY 100; 
PRED= (1.75* (10** (-~9)))+(4.87* (10**(-13))) *AGE; 


W= (AGE-2821.968) * (AGE-2821.968) ; 

SXX= 1928683* (202-1); 

SEPRED = SORT (6.204807* (10**(-19))) * SQRT(1+(1/202) +(W/SXX) ) ; 
U95= (PRED-1.645*SEPRED) ** (-1/3) ; 
L95=(PRED+1.645*SEPRED) ** (-1/3) ; 

AGE30=30; 

OUTPUT ; 

END ; 


PROC PRINT DATA= FIM4EXT;VAR AGE PRED L95 U95 ; 


PROC PLOT; 
PLOT L9O5*AGE='*' PRED*AGE=’ .’ U95*AGE='[’ AGE30*AGE=’.‘ /OVERLAY; 


DATA ONE( DROP=OBS); SET ALL1; 
FUME=EAGE/365; AGE=EAGE/365; 
PROC PRINT; VAR LOTID AGE FUME; 


DATA ONE; SET ALL1; IF AGE < 3 THEN DELETE; 


/* WITHIN-LOT MODEL */ 


PROC MIXED; 
CLASS LOTID; 
MODEL FUME = LOTID AGE*LOTID / NOINT S PREDICTED; 


/* FIXED EFFECTS REGRESSION MODEL */ 


PROC MIXED; 
CLASS LOTID; 
MODEL FUME = AGE / S PREDICTED; 
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/* RANDOM EFFECTS REGRESSION MODEL */ 


PROC MIXED ABSOLUTE; 
CLASS LOTID; 


MODEL FUME = AGE / §S PREDICTED; 
RANDOM INT AGE / CLG _ S SUB=LOTID TYPE=UN ; 
PROC PRINT; 


/* COMPUTE PRED SEPRED U95 L95 FOR RANDOM EFFECTS MODEL * / 


DATA TWO(DROP= OBS); SET ONE; BY LOTID; 
PRED=3 .81021477-0.14105394 * AGE; 
SEPRED=SQORT(( 0.070791) +(2*AGE* (-0.00450653) ) 
+ (AGE*AGE*0.00030661) +(0.05926819) ) ; 
U95=PRED+1.645*SEPRED; 
L95=PRED-1.645*SEPRED; 
PREDD=PRED*365; 
U95D=U95*365; 
L95D=L95*365; 
PROC PRINT; VAR LOTID AGE FUME PRED PREDD U95D L95D; 


/* GROUP SHELF-LIFE ESTIMATION * / 


DATA EXTEN7; 

DO AGE=(5000/365) TO (10000/356) BY (100/365); 
PRED=3 .81021477-0.14105394 * AGE; 
SEPRED=SORT(( 0.070791) +(2*AGE* (-0.00450653) ) 

+ (AGE*AGE*0.00030661) ) ; 
U95=PRED+1.645*SEPRED; 
L95=PRED-1.645*SEPRED; 

PREDD=PRED* 365; 
USS5D=U95*365; 
195D=L95*365; 
AGED=AGE* 365; 
OUTPUT ; 
END ; 
PROC PRINT; VAR AGE AGED PRED PREDD L95D : 


/* SHRINKAGE ESTIMATION * / 





DATA THREE; INPUT LOTID INTERCEP SLOPE; 
CARDS ; 

11142 -0.06196541 0.00479259 
11143 ~0.11611040 0.01007674 
11144 -0.16769671 0.01415533 
11145 -0.09516341 0.00494655 
11146 ~0.08786149 0.00599071 
11147 -0.21248452 Q0.01167817 
11148 -0.11397421 0.00853627 
11149 -0.07768438 -0.00123725 
41150 -0.15424392 0.00950736 
11151 -0.12749461 0.00861269 
11156 -0.15799910 0.00818934 
ALS -0.17588616 0.01418657 


oa 


11158 -0.13427054 0.00509910 
i159 ~0.23813672 0.01251881 
11160 -0.22564272 0.01415942 
11181 -0.16269825 0.00759863 
11182 -0.14882165 0,01126227 
Liss 0.17588730 -0.01074349 
11184 0.21035950 -0.01320753 
11185 0.44493213 -0.02624763 
11188 0.24783741 -0.01470002 
T1189 0.22284679 -0.01325530 
11190 0.18150475 -0.01253802 
Ni oe 0.22307183 -0.01355682 
LLL 2 0.20898367 -0.01170315 
11193 0.19615700 -0.01181337 
11194 0.17762001 -0.01152574 
L1195 0.16493381 -0.01078221 
20000 0.00000000 0.00000000 





PROC SORT; BY LOTID; 
PROC PRINT; 


DATA FOUR; SET THREE ; 
INTERCI=INTERCEP + 3.81021477; 
SLOPE1=SLOPE-0.14105394; 

PROC PRINT;VAR LOTID INTERC1 SLOPE1; 


DATA FIVE; SET THREE ; 
INTERC1=INTERCEP + 3.81021477; 
SLOPE1=SLOPE-0.14105394; 

DO AGE=(5000/365) TO (15000/356) BY (100/365) ; 
TSSL = INTERC1 + SLOPE1 * AGE ; BY LOTID; 
TSSLD=TSSL* 365; 

AGED= AGE*365; 


OUTPUT ; 
END ; 


PROC PRINT; VAR LOTID AGE AGED TSSL TSSLD 
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APPENDIX C FIGURES 


The figures shown in this appendix are described in Chapter 
II and IV are as follows 


Figure C.1 -Figure C.3 Standard Linear Regression Model I. 


A 
Figure C.4 -Figure C.6 Standard Linear Regression Model II. 
Figure C.7 -Figure C.8 Linearizing Transformation Model 1. 
Figure C.9 -Figure C.11 Linearizing Transformation Model II. 
Figure C.12-Figure C.14 Standard Regression Model Based on 

Fleet return data. 

Figure ‘C215 Weighted Least Square Model. 
Figure C.16-Figure C.17 Power Transformation Model. 
Figure C.18-Figure C.21 Individual fitted Models. 
Prgure “Cl22 Group shelf-life estimation. 


ye. 











FUME TIME (UNIT : DAY) 





0 2000 4000 6000 
PROPELLANT AGE (UNIT : DAY) 





Figure C.1. Fume time vs propellant age based on the 
standard linear regression model I. 
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Figure C.2. Group shelf-life estimation based on the 
standard linear regression model I. 
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Figure C.3. Residual analysis based on the standard linear 
regression model I. 
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Figure C.4. Fume time vs propellant age based on the 
standard linear regression model II. 
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Figure C.5. Residual analysis based on the standard linear 
regression model II. , 
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Figure C.6. Group shelf-life estimation based on the 
standard linear regression model II. 
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Figure C.7. Fume time vs propellant age based on the 
linearizing transformation model I. 
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Figure C.8. Residual analysis based on 
transformation model I. 
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Figure C.9. Fume time vs propellant age based on the 
linearizing transformation model II. 
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Figure C.10. Residual analysis based on the linearizing 
transformation model II. 
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Figure C.11. Group shelf-life estimation based on the 
linearizing transformation model II. 
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Figure C.12. Fume time vs propellant age based on the fleet 
return data. 
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Figure C.13. 
return data 
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Group shelf-life estimation based on the fleet 
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Figure C.14. Residual analysis based on the fleet return 
data 
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Figure C.15. Residual analysis based on the weighted least 
square model. 
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Figure C.16. Residual analysis based on the power 
transformation model. 
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Figure C.17. Fume time vs propellant age based on the power 
transformation model. 
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Figure C.18. Fitted within-individual models based on the 
standard linear regression model I. 
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Figure C.19. Fitted within-individual models based on the 
standard linear regression model II. 
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Figure C.20. Fitted within-individual models based on the 
fleet return data. 
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Figure C.21. Fitted within-individual models based on the 
linearizing transformation model I. 


74 


DAY) 


2 
Z. 
e 
3 
= 
7 
(2 
= 
= 
fe 


95e/0 ONE-SIDED LOWER PREDICTION LIMIT. 


PROPELLANT AGE (UNIT : DAY) 





Figure C.22. Group shelf-life estimation based on the random 
effects regression model. 
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